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ABSTRACT 

The early stages of dynamical evolution of planetary systems are often shaped by 
dissipative processes that drive orbital migration. In multi-planet systems, convergent 
amassing of orbits inevitably leads to encounters with rational period ratios, which 
may result in establishment of mean motion resonances. The success or failure of 
resonant capture yields exceedingly different subsequent evolutions, and thus plays a 
central role in determining the ensuing orbital architecture of planetary systems. In 
this work, we employ an integrable Hamiltonian formalism for hrst order planetary 
resonances that allows both secondary bodies to have hnite masses and eccentricities, 
and construct a comprehensive theory for resonant capture. Particularly, we derive 
conditions under which orbital evolution lies within the adiabatic regime, and provide 
a generalized criterion for guaranteed resonant locking as well as a procedure for 
calculating capture probabilities when capture is not certain. Subsequently, we utilize 
the developed analytical model to examine the evolution of Jupiter and Saturn within 
the protosolar nebula, and investigate the origins of the dominantly non-resonant 
orbital distribution of sub-Jovian extrasolar planets. Our calculations show that the 
commonly observed extrasolar orbital structure can be understood if planet pairs 
encounter mean motion commensurabilities on slightly eccentric (e ^ 0.02) orbits. 
Accordingly, we speculate that resonant capture among low-mass planets is typically 
rendered unsuccessful due to subtle axial asymmetries inherent to the global structure 
of protoplanet ary disks. 

Key words: planets and satellites: dynamical evolution and stability, celestial me¬ 
chanics, methods: analytical 


1 INTRODUCTION 

The orbital architectures of the Solar System’s planetary 
and satellite populations, as well as the currently known 
aggregate of extrasolar planets, exhibit numerous peculiar 
features. In principle, each specihc attribute of the galac¬ 
tic planetary census entails delicate constraints that inform 
the dominant mechanisms responsible for its inception and 
subsequent evolution. Therefore, theoretical interpretation 
of the observed state of planetary systems ultimately holds 
the key to understanding their origins. 

Among the most striking characteristics inherent to the 
orbital distribution of planets and satellites is the modest 
predisposition for orbital resonances, or loosely speaking, 
pairs of orbits with nearly rational periods. In the Solar 
System, the preference for commensurability among satellite 
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pairs (over randomly distributed orbits) is well-recognized, 
and was hrst pointed out more than half a century ago by 
Roy & Ovenden (1954) (see also Goldreich 1965 and Der- 
mott 1973). 

In the extrasolar realm, the radial-velocity sub-sample 
points to a mild over-abundance of resonant giant planets, 
typically at orbital radii in excess of ~ 1 AU (Wright et 
al. 2011; Winn & Fabrycky 2014). Meanwhile, the period- 
ratio distribution of an extensive collection of sub-Jovian 
planets discovered by Kepler and other surveys exhibits no¬ 
table enhancements at values slightly outside of exact 2:1 
and 3:2 commensurabilities (Fabrycky et al. 2014). Instances 
of multi-resonant chains, akin to the Galilean satellites, have 
also been found^. 

Gollectively, these observations imply that planet for- 


^ Examples of such systems iuclude GJ 876, Kepler-79 aud 
Kepler-223. 
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mation environments can be congenial towards emergence 
of resonant systems. However, the overwhelming majority 
of observed commensurabilities are unlikely to be truly pri¬ 
mordial in nature. Instead, resonances generally result from 
dissipative convergent evolution of pairs of orbits that fol¬ 
lows the initial phase of conglomeration (Goldreich 1965; 
Allan 1969, 1970; Sinclair 1970, 1972). 

There exists a multitude of physical mechanisms that 
may cause orbits to slowly approach one-another. Specih- 
cally, migration facilitated by tidal dissipation (Goldreich 
& Soter 1966; Peale 1976, 1986; Yoder & Peale 1981), in¬ 
teractions with a gaseous disk (Goldreich & Tremaine 1980; 
Grida, Sandor, & Kley 2008; Kley & Nelson 2012), and scat¬ 
tering of debris within a planetesimal swarm (Fernandez & 
Ip 1984; Malhotra 1993, 1995; Kirsh et al. 2009) are the 
most frequently quoted transport processes relevant to plan¬ 
ets and satellites. 

While convergent migration is required for resonances 
to congregate, this process alone does not guarantee success¬ 
ful resonant locking. In fact, the outcome of an encounter 
with a mean-motion commensurability depends on the in¬ 
trinsic parameters of the system, such as the planet-star and 
planet-planet mass ratios, the convergence rate, as well as 
the degree of pre-encounter orbital excitation. However, in 
spite of the uncertainty in the prospect of resonant capture 
itself, the post-encounter evolution of the system may de¬ 
pend critically on its result. 

As an example of resonances’ decisive nature, consider 
the fate of a pair of giant planets, at an epoch when the 
gaseous protoplanet ary disk has not yet dispersed. Although 
individually the planetary orbits would decay towards the 
central star along with the accretionary flow of the gas, suc¬ 
cessful capture of such objects into mean-motion resonance 
can reverse the overall inward drift and allow the planets 
to remain at large orbital radii^ (Masset & Snellgrove 2001; 
D’Angelo & Marzari 2012). As another example, assembly 
of sub-Jovian planets into resonances facilitates stabiliza¬ 
tion by the phase-protection mechanism (Greenberg 1977) 
and allows chains of planets to coherently migrate inwards 
without encountering each-other^ (Terquem & Papaloizou 
2007; Gresswell & Nelson 2008; see also Lee & Peale 2002). 
On the contrary, failure of objects to capture into resonance 
leads to an impulsive change in the system’s orbital proper¬ 
ties (Tittemore & Wisdom 1988, 1989, 1990) and may trigger 
large-scale dynamical instabilities (Tsiganis et al. 2005; Mor- 
bidelli et al. 2007; Batygin & Brown 2010). Thus, quantifi¬ 
cation of the capture probability of planets into mean motion 
resonance is central to understanding the early evolution of 
planetary systems. 

Over the last half century, numerous efforts have been 
made to develop a complete understanding of resonant cap¬ 
ture. Following the initial studies of Goldreich (1965); Gol¬ 
dreich & Soter (1966); Allan (1969, 1970); Sinclair (1972); 
Greenberg (1973), the hrst analytical formulation of the cap¬ 
ture probability within the framework of the circular re- 

^ Indeed, it is believed that this exact process is responsible for 
the retention of Jupiter and Saturn in the outer regions of the 
Solar System (Morbidelli & Crida 2007; Pierens et al. 2014). 

^ The long-term post-nebular stability of such systems is a sep¬ 
arate, non-trivial issue (Chambers, Wetherill, & Boss 1996; Ma- 
hajan & Wu 2014). 


stricted three-body problem (i.e. where a massive planet is 
assumed to orbit the central star on a circular trajectory 
and perturb a massless test-particle) was proposed^ by Yo¬ 
der (1973, 1979) and independently by Neishtadt (1975). 
Fully employing the use of adiabatic invariants, the formal¬ 
ism for analytic determination of capture probabilities was 
generalized by Henrard (1982) (see also Henrard & Lemaitre 
1983; Lemaitre 1984) and subsequently simplihed by Bor- 
deries & Goldreich (1984) (see also Peale 1986; Malhotra 
1988). Additional studies aimed at extending resonant cap¬ 
ture theory to the non-adiabatic regime have since also been 
undertaken by Friedland (2001); Quillen (2006); Ketchum, 
Adams, & Bloch (2011) and Ogihara & Kobayashi (2013). 

While the paradigm of the circular restricted three-body 
problem lends itself easily to theoretical analysis, it is not di¬ 
rectly applicable to the ultimately relevant issue of resonant 
capture of massive planets with non-circular orbits. Accord¬ 
ingly, to attack this more complicated problem, numerous 
authors have resorted to the use of numerical experiments. 
To this end, simulations that mimic the dissipative effects 
with hctitious forces (Lee & Peale 2002; Hahn & Malhotra 
2005; Terquem & Papaloizou 2007; Rein & Papaloizou 2010) 
as well as self-consistent hydrodynamical (Kley et al. 2005; 
Papaloizou & Szuszkiewicz 2005; Grida, Sandor, & Kley 
2008; Gresswell & Nelson 2008), self-gravitational (Moore, 
Quillen, & Edgar 2008; Moore & Quillen 2011) and tidal 
calculations (Mardling 2008; Lithwick & Wu 2012; Batygin 
Sz Morbidelli 2013a) have been performed. 

Although the published simulations represent impor¬ 
tant advancements in the understanding of planet formation 
and evolution, such calculations continue to leave analytical 
development to be desired for two reasons. First, numerical 
experiments are typically tailored towards particular sys¬ 
tems with specihc choices of physical parameters. Without 
the knowledge of underlying theoretical relationships, it be¬ 
comes difficult to translate the obtained results to another 
case without additional modeling. Second, despite substan¬ 
tial advances in computational technologies, high-resolution 
numerical experiments remain much too computationally 
expensive to practically serve as a replacement for theory. 

Gonsequently, in order to explain why planetary sys¬ 
tems possess the orbital architecture that they have, it is 
necessary to generalize the existing resonance capture the¬ 
ory to the physical domain of the unrestricted elliptic three- 
body problem. This is the primary purpose of this work. Ul¬ 
timately, with a more complete model in hand, we shall de¬ 
lineate parameter regimes within which resonant capture is 
expected to be successful, and aim to understand why res¬ 
onant orbits are neither ubiquitous nor absent within the 
galactic planetary census. 

The paper is structured as follows. In section 2, we gen¬ 
eralize the existing hrst order resonance capture theory to 
allow for both planets to have masses as well as hnite eccen¬ 
tricities. Additionally, we obtain a criterion which dictates 
whether or not the dissipative evolution lies in the adiabatic 
regime. In section 3, we examine some immediate conse¬ 
quences of the model. Particularly, we construct maps of 


^ It should be noted that the concept of probabilistic capture 
within the framework of spin-orbit resonances was hrst introduced 
by Goldreich & Peale (1966). 
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Figure 1. Equilibria of Hamiltonian (13). The three roots of 
equation (16) are shown as functions of the resonance proxim¬ 
ity parameter. The elliptic equilibrium point that is always real 
and negative is shown in blue. The elliptic and hyperbolic fixed 
points that only exist above the bifurcation value of (5 ^ (5c are 
shown in black and red respectively. Extrinsically driven conver¬ 
gent evolution of the orbits corresponds to a gradual increase in 
(5. 


capture probabilities relevant to adiabatic encounters and 
derive specific adiabatic thresholds for disk-driven migra¬ 
tion, planetesimal driven migration, and tidal evolution. 
Subsequently, we analytically consider the orbital evolution 
of Jupiter and Saturn while submerged within the proto¬ 
solar nebula, and derive critical eccentricities above which 
sub-Jovian extrasolar planets fail to capture into mean mo¬ 
tion resonances. We summarize and discuss our results in 
section 4. 


2 ANALYTICAL THEORY 

The calculation we aim to perform is perturbative in na¬ 
ture. Accordingly, we begin with a basic description of the 
dynamics. It is standard practice in celestial mechanics to 
treat planet-planet interactions as perturbations to Keple- 
rian motion. As such, planet-planet potential is typically ex¬ 
panded as a Fourier series in the orbital angles and a power 
series in the planetary eccentricities and inclinations (Laskar 
Sz Robutel 1995; Ellis Sz Murray 2000; Laskar & Boue 2010). 

In the vicinity of a given resonance, it is sensible to av¬ 
erage over all short-periodic terms (i.e. those that vary on an 
orbital timescale) in the expansion, and retain only slowly 
varying (compared to the orbital period) harmonics, which 
include those associated with the resonant interaction (Mur¬ 
ray Sz Dermott 1999). This procedure yields a simplified 
Hamiltonian that approximates the full Hamiltonian near a 
commensurability. 


2.1 Model Hamiltonian 

In its most rudimentary form (i.e. to first order in planetary 
masses and eccentricities), the planar k \ k — 1 resonant 
Hamiltonian reads: 

^ _ QM^mi QM^m2 Qmim2 

2ai 2a2 a2 
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fres^i cos(A:A2 - {k - l)Ai - mi) 


+ /reie2 cos(A;A2 - (k - l)Ai - 1172) 


( 1 ) 


In the above expression, k is an integer greater than unity, Q 
is the gravitational constant, is the mass of the central 
star, m is the planetary mass, a is the semi-major axis, e 
is the eccentricity, A is the mean longitude, m is the lon¬ 
gitude of perihelion, and /res’s are functions of order unity 
that (weakly) depend on the semi-major axis ratio (ai/a 2 ) 
(see Callegari Sz Yokoyama 2007 for explicit expressions). 
The subscripts 1 and 2 denote the inner and outer planets 
respectively. 

Upon inspection it is immediately clear that the Hamil¬ 
tonian (1) contains two resonant arguments that appear in 
two separate harmonics. Although at first glance this Hamil¬ 
tonian appears to comprise a typical two degree of free¬ 
dom system, it was shown to be integrable by Sessin (1981); 
Sessin & Ferraz-Mello (1984). The published calculations of 
Sessin & Ferraz-Mello (1984) are rather involved, as they 
were performed within the framework of the Hori (1966) 
perturbation method. However, the demonstration of inte- 
grability was greatly simplified by Wisdom (1986) and Hen- 
rard et al. (1986), who showed that instead of constructing 
formal perturbation series, it is possible to consolidate the 
two harmonics into a single term with the aid of a reducing 
canonical transformation^ that corresponds to a rotation in 
phase-space (see also Henrard & Lemaitre 2005). 

By employing the reducing transformation, we shall 
cast the Hamiltonian into a form that is synonymous to 
a first-order Andoyer Hamiltonian, also known as the sec¬ 
ond fundamental model for resonance (Henrard & Lemaitre 
1983; Morbidelli 2002). In this work, we will not spell out 
the derivation in excruciating detail, as it is available else¬ 
where (Sessin & Ferraz-Mello 1984; Wisdom 1986; Henrard 
et al. 1986; Ferraz-Mello 2007; Batygin & Morbidelli 2013b; 
Deck, Payne, Sz Holman 2013; Delisle, Laskar, Sz Correia 
2014). Instead, we shall restrict ourselves to sketching out 
the important steps. In both notation and substance, we 
shall closely follow the derivation outlined by Batygin & 
Morbidelli (2013b). 

As a first step, we define the canonical variables 

K: = Al + k = Al 

rv 

0 = A 2 /A: 0 = kX2-{k- l)Ai 

Xi = eiy/[K[iCOs{-mi) yi = ei^/[K\iS\n{-Wi) (2) 

where A^ = rrii^/QM^ai and the square brackets [ ] de¬ 
note an evaluation at nominal resonant semi-major axis (i.e. 
[a ]2 = (k/(k — l))^'^^[a]i). Upon expanding the Keplerian 
part of the Hamiltonian to second order in A around [A] 
(which corresponds to a perturbation that is first order in 
e), neglecting semi-major axis variation in the disturbing 
part of the Hamiltonian (since it is already first order in e), 
and dropping constant terms, the expression (1) takes the 
form: 

H = 3[h]iik - 1)/C0 - 5 {[h]iik - if + [h]2C)e^ 


^ Evidently, a transformation of this sort was first proposed by 
Poincare (1899). 
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Figure 2. Evolution of the phase-space portrait of the Hamiltonian during convergent orbital migration. Panels A, B and C depict 
level curves of equation (13) for 5 = —1, 0, and 1 respectively. The equilibrium points are shown with small dots and are color-coded 
in the same way as in Figure (1). While libration of the critical angle is possible in panel A, formal resonant trajectories only exist in 
panels B and C, where the separatrix is shown as a thick red curve. Across the panels, three values of the phase space area (adiabatic 
invariant) are accentuated: sub-critical (labeled J' and shown as a shaded blue region), critical (labeled J7c and emphasized with a green 
shade in panel B), and super-critical (corresponding to the area occupied by the outer branch of the separatrix on panel C). Trajectories 
encircling the critical and super-critical areas are shown on panel A as thickened red lines. A circulating orbit engulfing a super-critical 
area is similarly shown in panel B. While smooth adiabatic entry into resonance is possible for initially circulating trajectories occupying 
a phase-space area equal to, or smaller than J7c, those occupying a larger area must cross the separatrix to be captured, as shown on 
panel C. On panel B, the minimal and maximal excursions of the action along the separatrix (which correspond to the resonance width 
at the onset of resonance) are labled and p^^ respectively. 


- {axi + ^X 2 ) cos( 6 >) + {ayi + ^y 2 ) sin( 6 >). (3) 

Here, [h] = N/[A] = l/(m[a]^) is the inverse moment of 
inertia of a circular orbit, and 


_ f^es 

"" ~ [A]i y[Al7 

... 

^ [All ^[A]^ ^ ^ 

are constants that encapsulate the strengths of the individ¬ 
ual harmonics present in the Hamiltonian. 

Hamiltonian (3) is independent of the angle k, mean¬ 
ing that the action /C is a constant of motion. However, the 
number of degrees of freedom is still too great for integra- 
bility. Defining a canonical rotation (Wisdom 1986; Henrard 
et al. 1986; Ferraz-Mello 2007) 


ui = 


U2 = 


axi + 13x2 
•JcP- + /32 
I3xi — ax2 


vi = 


V2 = 


ayi + Py2 
+ /32 
Pyi - ay2 


s/oF+W -Jo? + /32 

and the associated action-angle coordinates 


(5) 


This reduction identifies a second integral of motion, that is 
^ 2 . 

Employing a contact transformation, we define a final 
set of action-angle variables 

'0 = 01 + 6 > 

Q = ^ cr = 6 > ( 8 ) 

which renders H. independent of cr. Accordingly, identifying 
Gl as the final conserved quantity and dropping constant 
terms, the Hamiltonian is cast into an integrable form: 

u = -‘i {[h]2pn - [h]i{k -1) (/c - (fc - i)n)) ^ 

-l{lh]i{k-lf + [h]2p)'i>^ 

— cos('0). (9) 

The constants that precede each of the terms in the 
Hamiltonian are not truly independent. Thus, as a conclud¬ 
ing step in the preparation of the Hamiltonian, we shall scale 
the action and time, in order to introduce a single resonance 
proximity parameter^ 5. To maintain symplecticity, we scale 
all of the actions present in the Hamiltonian, as well as the 
Hamiltonian itself by the same constant factor, 77 : 


Ui — \/2^z cos(0i) V — \/2$z sin(00, ( 6 ) 


ii = UIp ^ = T^/77 iC^ICIr] Q = Q/ry. (10) 


the Hamiltonian is reduced to a form that contains only a 
single harmonic: 

U = 3[h]iik - 1)/C0 - 5 {[h]iik - if + [h]2p)e^ 

- ya2 +^2y2^cos(^i +6>). (7) 


We choose the expression for 77 such that the coefficient up¬ 
front the term that is quadratic in action is half® as big as 


® This is a slightly different scaling from that employed by Baty¬ 
gin & Morbidelli (2013b). 
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Figure 3. Area index and capture probability. Panel A shows the phase-space area occupied by the inner and outer branches of the 
separatrix as a function of the proximity parameter. In the adiabatic regime, the initial phase-space area occupied by a circulating 
trajectory far away from resonance matches onto the area of the outer branch of the separatrix, evaluated at the proximity parameter, 
(^enc, at which the resonant encounter occurs. Panel B shows the resonance capture probability as a function of the initial action. Note 
that the capture probability diminishes rapidly for values of the initial action that exceed the critical value below which capture is certain. 


the one that precedes the harmonic: 
r? = ■ 


9{[h]i{k - 1)2 + [h]2k^)^ 

(/re^i)^mi + {f^el)‘^{k/{k - 
9M2 {k^mi {k — 1 ) 2 / 3 /^ 4 / 37722 )^ 


= 17111712 \/GM^[a]2 

1 1/3 


( 11 ) 


Subsequently, we change the unit of time to 3r]{[h]i{k — l)‘^-\- 
[h\2k‘^)/2, and divide H by the same factor. Writing the 
proximity parameter as 

(5 = -{[h]2k‘^{3 + 2Q) + [h]i{k - 1) 

X (3A; - 2/C + 2(/c - 1)Q - 3)) 

X {[h]iik-lf+ [h]2C)-\ (12) 


the Hamiltonian takes on the familiar form: 

= 3(5 + 1)^ - - 2^/2!'cos(V>). 


(13) 


With the exception of the signs^, this is the Hamiltonian 
considered by Yoder (1973, 1979); Neishtadt (1975); Hen- 
rard (1982) and Borderies & Goldreich (1984). Note that 
this expression for the Hamiltonian is written in terms of 
dimensionless variables. Thus, the quantity 77 (which has 
units of angular momentum) encapsulates all of the infor¬ 
mation regarding how the dynamics scales with mass ratios 
and physical sizes of the orbits. 


2.2 Equilibria of the Hamiltonian 

Some insight into the properties of the Hamiltonian can be 
obtained by considering its equilibria. Defining the cartesian 
coordinates 

p = cos( 7 /i) q = \/2^ sm{ip), (14) 


^ Signs only determine if the resonant equilibrium point resides 
at 7 /; = 0 or 7 /; = TT, and are thus unimportant to the discussion 
at hand. 


the Hamiltonian can be written as: 


n = 3{s + i) 






2p. 


(15) 


It can be understood from equation (15) that the equilib¬ 
rium points must reside on the p—axis, since 7 /; = 0, tt for 
d'H.jdil) — 0. Thus, setting ^ = 0, the equilibrium equation 
reads: 


3{S + l)p — — 2 = 0. (16) 

The cubic equation (16) admits three roots, one of 
which is always real and negative, while the other two are 
real (and positive) only above a critical bifurcation value of 
5, which can be computed by setting the two positive roots 
equal to each other. For Hamiltonian (13), the bifurcation 
occurs at (5c = 0 . Strictly speaking, S ^ Sc defines a condi¬ 
tion for the existence of resonance, since a homoclinic curve 
(i.e. separatrix), necessary for identihcation of resonant tra¬ 
jectories, only exists for values of S greater than or equal 
to Sc (see e.g. Delisle et al. 2012 ). To this end, it is worth 
noting that although the functional form for S is compli¬ 
cated, it monotonically grows as orbits approach each-other 
convergently. Thus, an increase in the value of S unilaterally 
signals evolution into the resonance. 

The three roots of equation (16) are shown as functions 
of S in Figure (1) and can be classihed as follows. The neg¬ 
ative root (shown with a blue line) always corresponds to 
a stable (elliptical) equilibrium point and for (5 ^ (5c, is en¬ 
veloped by resonant trajectories. The positive root, whose 
value is closer to zero (shown with a black line), also cor¬ 
responds to a stable hxed point but lies at the center of 
the inner circulation region of the phase space portrait. The 
larger positive root (shown with a red line) corresponds to 
an unstable (hyperbolic) equilibrium. For reference, phase- 
space portraits of the Hamiltonian for (5 = — 1, (5 = (5c = 0 
and (5=1 are depicted in Figure (2), where the hxed points 
are color-coded in the same way as in Figure (1). 

Note that the unstable equilibrium dehnes a point 
where the inner and outer branches of the separatix join. 
For this reason, its value is particularly important for the 
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Figure 4. The dimensionless libration frequency associated with 
the resonant equilibrium point. The frequency is minimized at a 
value of the proximity parameter slightly below 6c and increases 
slowly and monotonically for 5 > (5c. In this work, we adopt the 
value of cu at (5 = (5c as the characteristic value. 


evaluation of capture probabilities. Written explicitly, the 
corresponding expression for the unstable fixed point reads: 

Peq — — ty/s ^\/— 6(^6(^6 + 3) + 3) + 

X ^2 ^ ^/^6(S{S + 3) + 3) + , (17) 

where i — \/^. Recall that this expression is purely real 
only for (5^0. 

2.3 Conditions for Guaranteed Capture 

As already stated above, within the context of Hamiltonian 
(13), the slow convergence of two initially non-resonant or¬ 
bits towards a mean motion commensurability can be en¬ 
visioned as a gradual increase in 6 from values below 6c 
to values above 6c (Peale 1986). Provided that dissipative 
processes responsible for secular changes in planetary archi¬ 
tecture act on much longer timescales than the dynamical 
timescales associated with resonant motion, we can define 
an adiabatic invariant (Henrard 1982; Neishtadt 1984) 

J = ^ ^ d'lp = ^ p dq, (18) 

which is conserved as long as the trajectory does not 
encounter a homoclinic curve (Lichtenberg & Lieberman 
1983). Physically, J corresponds to the area occupied by 
the orbit in phase space. 

For a given value of (5, the conservation of J dictates the 
energy level on which the trajectory resides. This is demon¬ 
strated in Figure (2) where a set of orbits characterized by 
the same J are shown on phase-space portraits with differ¬ 
ent values of 6. The condition for guaranteed capture into 
resonance stems directly from the fact that the adiabatic in¬ 
variant is preserved as long as the orbit does not encounter 
a separatrix. 

Of the three thickened circulating trajectories shown in 
panel A of Figure (2), consider the one with the smallest 


radius (i.e. the orbit whose engulfed area is shaded blue and 
labeled J7). This orbit is bound to become a resonant tra¬ 
jectory, (specifically the one corresponding to J in panels B 
and C), because the phase-space area it occupies is smaller 
than the area occupied by the separatrix at (5 = (5c, (we 
shall call this area J'c), and the total phase-space area occu¬ 
pied by the resonance is minimized when the separatrix first 
appears (Henrard 1982; Morbidelli 2002). Furthermore, be¬ 
cause the area occupied by resonant trajectories only grows 
with increasing (5, the libration amplitude of the trajectory 
in question shall only decrease with time (Yoder 1979; Peale 
1986). For reference, J7c is shown in green shade in Figure 
(2B). 

Let us calculate J7c. Setting 6 = 6c, from equation (17) 
we obtain p*^ = 1. Accordingly, the value of the Hamiltonian 
that corresponds to the separatrix, Tic, is given by 

H£ = 3(<5c + 1)-^--^-2p:, = -|. (19) 


With a value of l-Lc in hand, the relationship between ^ 
and '0 along the separatrix at (5 = (5c can be obtained® from 
equation (13). 

Integration along the homoclinic curve shows that the 
phase-space area occupied by the separatrix is given by 


Jc = 6 


arcsin 


1 

(Peq)^''^ 



= 67r. 


( 20 ) 


The same procedure yields the area occupied by any tra¬ 
jectory (identified by the value of 7i) for any value of 6. 
The evaluation of is particularly trivial in the limit where 
(5 (5c. In this case, the phase-space portrait is composed 

of concentric circles, centered on the stable fixed point near 
the origin. Therefore, provided some value of ^ = p^/2, we 
have 


= Tvp^ = 27r 


( 21 ) 


Naturally, this limit is quite relevant for convergently mi¬ 
grating planets or satellites, as it represents the dynamical 
state of the system far away from mean motion commensu¬ 
rability. 

The above discussion indicates that the evaluation of ^ 
far away from resonance is sufficient to gauge whether or not 
guaranteed resonant capture will take place. The calculation 
of ^ itself requires the specification of planetary eccentric¬ 
ities, longitudes of perihelia as well as planet-star mass ra¬ 
tios. Provided the appropriate quantities, we can write the 
expression for ^ in terms of orbital elements and thus for¬ 
mulate a practically useful criterion for guaranteed adiabatic 
capture within the framework of the unrestricted three-body 
problem. Specifically, we have: 



where ( = {{k — 1)/A;)^/®. 


® Alternatively, an equivalent relationship between p and q can 
be obtained from equation (15). 
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Figure 5. Resonant libration period as a function of planet-planet and secondary-primary mass ratios. Panel A shows the dependence 
of the libration period on how mass is distributed among the two secondary bodies. Panel B depicts the libration period as a function 
of the total secondary mass to the central mass. Clearly, rnb is independent of the planet-planet mass ratio to a good approximation. 
Relationships corresponding to the 2:1, 3:2 and 4:3 resonances are shown as red, blue, and green lines respectively. 


The quoted criterion requires knowledge of the plane¬ 
tary eccentricities as well longitudes of perihelia, which is 
rather inconvenient since these quantities generally depend 
on the secular dynamics of the planets outside of the reso¬ 
nant domain. Moreover, neither T nor $2 (identified in equa¬ 
tion 5) are individually conserved in the secular domain, and 
large scale variations of the two quantities are indeed pos¬ 
sible in principle (see Batygin & Morbidelli 2013b). On the 
other hand, the angular momentum deficit 

A = r]-^ + <i>2 = mi CQM ,Hi y + m 2 y/QM ,[aja y, (23) 

which is more readily interpreted, is conserved on the secular 
domain. 

Because $2 ^ 0, we have: T ^ A/r]. Taking advantage 
of this, we can formulate an excessively strong, but very 
simple version of the criterion for guaranteed capture: 

- ^ 3 . (24) 

T] 

It is important to note that the maximal extent of orbital 
excitation that allows for guaranteed capture is a function 
of the secondary to primary mass ratios. In particular, equa¬ 
tions (22) and (24) imply that the maximum allowed eccen¬ 
tricities decrease with decreasing mass-ratio, a relationship 
that may be anticipated intuitively. 


2.3.1 Circular Restricted Problem as a Special Case 


Although expression ( 22 ) is derived here for the first time, 
a related criterion for guaranteed adiabatic capture is well 
known within the framework of the circular restricted three- 
body problem (see Murray & Dermott 1999). Provided that 
equation ( 22 ) is a more general criterion, it should reduce to 
the restricted limit as one of the planetary masses is taken 
to vanish, while the orbit of the other planet is assumed to 
approach a circle. 

Considering interior resonances first, we set mi = 0 and 
62 = 0 in equation ( 22 ) to obtain: 


61 ^ 


m2 frel 

3M* /c2/3(/._ 1)4/3 


(25) 


Similarly, for the case of exterior resonances we set m 2 = 0 
and 61 = 0 , which yields 


62 ^ Vd 


mi /i 


( 2 ) 

res 


1 1/3 


3M* 


(26) 


Both of these expressions are in direct agreement with those 
derived by assuming a circular restricted formalism from the 
beginning (Peale 1986; Malhotra 1988). 

It is noteworthy that within the framework of the re¬ 
stricted problem, the aforementioned capture criteria for in¬ 
terior and exterior resonances are obtained from different 
Hamiltonians. Specifically, the analysis typically begins with 
a variant of equation (1) where only a single resonant har¬ 
monic is identified. Subsequently, the corresponding expres¬ 
sions are individually cast into the form of Hamiltonian (13) 
and a calculation of capture probabilities is carried through. 
In this work, we have utilized the reducing transformation 
(5) to derive the previously known criteria as limiting cases 
of a more comprehensive formalism. Our analysis thus gen¬ 
eralizes previous results. 


2.4 Probabilistic Capture 

If the criterion (22) is not satisfied, resonant capture is not 
certain. This is because at sufficiently high values of the 
action T, the phase-space area occupied by the orbit exceeds 
the critical area occupied by the separatrix at the inception 
of the resonance i.e. 2T > Jc at (5 = (5c. Consequently, unlike 
the “smooth” transition to resonance discussed above, in 
this case the trajectory must cross the separatrix in order 
to reach the resonant domain. However, passage across a 
critical curve is a fundamentally probabilistic process, and 
may advect the trajectory into the inner circulation region 
of the phase-space portrait. 

To determine the probabilistic outcome of an encounter 
of a > jTc trajectory with a homoclinic curve, one must 
first relate the non-resonant initial conditions to the value 
of the proximity parameter at which the transition will take 
place. Utilizing the conservation of the adiabatic invariant 
(equation 21), we can connect the value of the action To far 
away from resonance (i.e. at (5 ^ (5c, where the trajectory in 
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phase-space can be approximated by a circle) to the value 
of the proximity parameter at which the encounter will take 
place Sene- This can be done by matching the phase space 
area occupied by the trajectory to that engulfed by the sep- 
aratrix at Sene- The appropriate expression reads (Henrard 
1993): 




3CMi,(2(k — 1)^1712 + 


1 2/3 


( {frea ) ^^2 + C (/rei) ^”^1) 

X ( (/re/)^ei + (/rei)^e2 + 6162 COs{Aw) 

1 


— ^ ^((^enc "h 1) 

V(PSq)^ - 1 




TT 

+ 2 


(27) 


Recalling that Peq is a single-valued function of 5 (see equa¬ 
tion 17), the above expression unequivocally links the start¬ 
ing condition of the system to the phase-space portrait at 
which the encounter takes place. Figure (3A) shows Jj 21 ^ — 
To as a function of (5enc, obtained from equation (27). 

With the above relationship at hand, we can evaluate 
the probability of resonant capture for arbitrary initial con¬ 
ditions. Generally, adiabatic capture probability, P, is given 
by the ratio of the rates of changes of the phase-space areas 
occupied by the resonant domain of the phase-space por¬ 
trait, to that occupied by the outer branch of the separatrix 
(see Figure 2C). Following Borderies h Goldreich (1984), the 
explicit expression® for V that corresponds to the Hamilto¬ 
nian (13) reads: 


(28) 


Figure (3B) depicts the resonant capture probability as 
a function of Tq. It should be noted that in the case of 
probabilistic capture, the transition through the separatrix 
necessarily breaks the conservation of J. This means that 
the resonant encounter will lead to an impulsive change in 
orbital parameters, irrespective of whether capture is suc¬ 
cessful. To this end, note that passage through the resonance 
(see Figure 2G) diminishes the overall phase-space area, im¬ 
plying a less dynamically excited post-encounter state. 


V = 2 

11 ^ 

-1 

2arcsin ((plq)-^''^) 



2.5 The Adiabatic Threshold 

The entirety of the above discussion assumes that the dissi¬ 
pative evolution of the proximity parameter is so slow that at 
any point, a “frozen” (i.e. constant (5) representation of the 
dynamics provides a good approximation to the real evo¬ 
lution (Wisdom 1985; Henrard & Garanicolas 1990). This 
approximation holds in a regime where the resonant libra- 
tion frequency greatly exceeds the extrinsic resonance cross¬ 
ing rate. Within the framework of the adopted model, both 
quantities can be analytically derived, thereby providing an 
explicit criterion for adiabatic dynamics. 


® It should be understood that within the framework of this ex¬ 
pression, a value of V that exceeds unity simply corresponds to 
certain capture. 


To calculate the resonant libration frequency, we ex¬ 
pand the Hamiltonian (13) to second order in T and ^ 
around the equilibrium point ([T],7r) that corresponds to 
the definite negative root of equation (16). Dropping con¬ 
stant terms, we have: 


+ ( 3 + 35-2m+ 7^) ('I-[#]). (29) 


The fixed-point condition dictates that in the above expres¬ 
sion, the term linear in T must vanish. This yields a relation¬ 
ship between the proximity parameter and the equilibrium 
action: 


5 = 



(30) 


Introducing new canonical variables (e.g. Batygin h 
Morbidelli 2013b): 


T = (T- [T]) 


'0 = ('0 — Tt) 


/y 2 + 4 [^y^" 

V ) 

4 v ^[^]2 J 


(31) 


we obtain a Hamiltonian that is synonymous to that of a 
simple harmonic oscillator: 


n = + (32) 

where 


CJ = 


2 + 4^/2[T]3/2 


[T] 


(33) 


is identified as the resonant libration frequency (Lichtenberg 
& Lieberman 1983). 

Utilizing equation (30), cu can be expressed as a func¬ 
tion of the proximity parameter, S, and this relationship is 
shown in Figure (4). It is trivial to show that the value of 
cj is minimized at (5 = — 1 (corresponding to [T] = 
and quantitatively evaluates to cjmin = 2^/^\/3. For larger 
values of (5, cj slowly increases. For example, at the onset of 
resonance (i.e. at 5 = 0 when the separatrix first appears), 
cj = 3. Recalling the scalings introduced in equations (10) 
and (11), we adopt the value of the libration frequency at Sc 
to formulate a practically useful criterion for adiabatic evo¬ 
lution^®. Specifically, in terms of physical parameters, the 
corresponding resonant libration period is given by: 


An {k — 1)^/® 

3\^ 3V3A;5/9 

/_M*_\ 

\{{k - l)/fc)l/®(/resU”T’l + (/re’'i)^m2/ 


Non-adiabatic evolution will cause the trajectory to jump to 
the inner circulation region without being trapped into resonance. 
So, it is sensible to consider the libration frequency at a value of 
6 at which the inner circulation region first appears i.e. Sc- At 
this value of (5, a; exceeds its minimum value by ~ 10%. 
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Figure 6. The adiabatic threshold as a function of plane-planet and planet-star mass ratios. As in Figure (5), panel A depicts the 
dependence of the adiabatic criterion on the distribution of mass among the planets, while panel B shows its —4/3 power-law dependence 
on the ratio of the total planet mass to the central mass. In addition to the theoretical relationships, numerically computed adiabatic 
threshold data corresponding to equal-mass bodies in a 2:1 resonance (Ogihara & Kobayashi 2013) are shown with blue points. Curves 
corresponding to the 2:1, 3:2 and 4:3 resonances are shown as red, blue, and green lines respectively. 




1/3 


^ {{k — l)k‘^y/^mi + (/c — l)m2 J ' 

It is instructive to examine the dependence of rub on 
the distribution of mass between mi and m 2 . The relation¬ 
ship can be evaluated by setting mi = /i(l — M* and 
Tn 2 = where fi — {mi + 7712 )/M*, and examining 

the dependence of rub on Figure (5 A) shows rub, appro¬ 
priately scaled by the mass ratio parameter and the outer 
planet’s orbital period for 2:1, 3:2, and 4:3 resonances. As 
can be immediately gleaned, the libration period is indepen¬ 
dent of ^ to an excellent approximation^^. Thus, we can sim¬ 
plify the above expression somewhat by replacing m 2 with 
the cumulative mass of the planets and setting mi = 0: 


' 7 ”lib — 


An 

( ) 

2/3 

1 

■(3(AsT)-^/^' 

3 [n ]2 ' 

7771 + 7772 

(fc5(fc- l))l/9 


(35) 


Figure (5B) shows the resonant libration period as a function 
of the mass-ratio. 

In order to obtain a condition for adiabatic evolution, 
we require an additional piece of information: the character¬ 
istic timescale on which dissipative forces will carry the sys¬ 
tem across the resonance width, in absence of planet-planet 
interactions (Friedland 2001). As a first step, we must de¬ 
fine the resonance width, AT. The maximal excursion in T 
at (5 = 0 corresponds to the difference between the points 
at which the critical curve crosses the p—axis. Noting that 
one of the crossing points is the bifurcated equilibrium, from 
equations (17) and (19), we have: 


AT = 


{py^ {pUf 


r] = dry. 


(36) 


2 2 

where p* = — 3 signifies the negative intersection point of 
the separatrix. 

From conservation of the integrals 1C and Cl (see equa¬ 
tions 2 and 8), it follows that the maximal resonant variation 
of the first Poincare momenta is 

AAi = -(/c- 1) ATry 


Specifically, the fractional variation is of order ~ 1%. 


AA 2 = k AT ry. 


(37) 


The quantity of interest for our calculation is the ratio of 
mean-motions: 

^ ^ ^ A 3 (Ax - [A]i) _ 3 (A2 - [A]2) \ ^ 33) 

ni [n]i V [A]i [A]2 J 

Accordingly, the above expressions define a resonance width 
in x: 


n]i V ' ^ [A]i 

k ([A]i + [A]2) — [A]2 


= —1277 


[A]i [A] 2 


(39) 


Let us now turn to dissipative evolution. A conventional 
way to parameterize extrinsically facilitated drifts in semi¬ 
major axes is to adopt the following form (Lee & Peale 2002): 


dai _ ai 
dt Ti 


da2 _ a2 
dt T2 


(40) 


From these expressions, it follows that the associated rate 
of change in the ratio of mean motions is: 


dx 

dt 


3 (ti - T 2 ) 772 

2 T1T2 ni ’ 


(41) 


By replacing the derivative on the LHS by a fraction of finite 
differences, and equating the change in x to the resonance 
width, we obtain a corresponding segment of time required 
for dissipative forces to carry the orbits across the resonance 
(Friedland 2001; Goldreich & Schlichting 2014): 

^ 877 ^ ~ . (42) 


T 1 T 2 


[A]i [A] 2 


The condition for adiabatic evolution is thus defined: 
the characteristic timescale of orbital convergence. At, must 
exceed the libration period rub. In terms of physical param- 
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ei 


Figure 7. Conditions for guaranteed adiabatic capture of Jupiter 
and Saturn into the 2:1 mean motion resonance. Each of the plot¬ 
ted curves corresponds to a distinct assumed difference in the 
apsidal lines, Atu. The ei — 62 parameter space associated with 
certain capture is clearly minimized at Atu = tt and maximized 
at Atu = 0. 


eters, the adiabatic criterion reads: 



Similarly to the preceding discussion of the expression 
for Tiib, equation (43) is approximately independent of how 
mass is distributed between the secondary bodies. Instead, 
the criterion depends sensitively on the ratio of the total 
planetary mass to the mass of the central object. Accord¬ 
ingly, in the same spirit as equation (35), we can obtain a 
simplified expression for the adiabatic criterion: 


Tlib ^ TT (n — T2) 

At 2 [n]2 T 1 T 2 




mi + m2 
< 1 . 


4/3 


(44) 


We note however, that dependence on the planet-planet 
mass ratio can still enter equations (43) and (44) through 
particular expressions for n and T 2 . 

The analytical expression for the adiabatic threshold 
derived above exhibits the characteristic —4/3 power-law de¬ 
pendence on the mass ratio. Unsurprisingly, the same power- 
law holds for the restricted problem (Friedland 2001; Quillen 
2006; Goldreich & Schlichting 2014), and has been observed 
to arise in numerical N-body simulations (Ketchum, Adams, 
& Bloch 2011; Ogihara & Kobayashi 2013). Figure ( 6 A) 


shows the dependence of the adiabatic criterion on the distri¬ 
bution of masses, while the dependence on the cumulative 
secondary to primary mass ratio is shown in Figure ( 6 B). 
Also over-plotted on panel B are numerical estimates of the 
adiabatic threshold for two equal-mass bodies correspond¬ 
ing to the 2:1 resonance obtained by Ogihara & Kobayashi 
(2013). Clearly, the agreement between theory and simula¬ 
tion is more than satisfactory. 


3 RESULTS 

With the crucial features of the analytical theory defined, we 
are now in a position to examine the generic consequences 
of the model. Additionally, we shall aim to understand how 
observations of exoplanetary systems and simulations of the 
early dynamical evolution of the Solar System fit into the 
framework of the developed formalism. We shall begin by 
outlining some generalities. 


3.1 General Results 

3.1.1 Capture Probability Maps 

The formulae presented in the previous section allow for the 
evaluation of P at a computationally negligible cost. Tak¬ 
ing advantage of this, we are provided with an opportunity 
to rapidly delineate the probability of capture in a given 
section of parameter space, a task that would be impossi¬ 
ble to accomplish with brute-force numerical simulations. To 
this end, the only quantity we need to specify is ^ far away 
from resonance. However, ^ itself represents a rather compli¬ 
cated combination of orbital and physical parameters, and 
involves the specihcation of planet-star and planet-planet 
mass-ratios, orbital eccentricities, as well as the difference 
in the longitudes of perihelia. 

Although the specihcation of mass-ratios and eccentrici¬ 
ties is unavoidable, one would ideally like to be agnostic with 
respect to specifying the difference in longitudes of perihelia 
of the orbits, Avj. Unfortunately, the outcome of resonant 
encounters is not insensitive to Avj. As an example, consider 
a pair of planets with star-planet mass ratios equal to that 
of Jupiter and Saturn in the vicinity of a 2:1 commensura- 
bility (this example will also be relevant to the discussion 
presented below). Utilizing equation (22), we can delineate 
curves that correspond to guaranteed capture on a ei vs 62 
digram, as shown in Figure (7). Clearly, parameter space as¬ 
sociated with certain capture is minimized for Avj — tt and 
maximized for Atu = 0. 

Thus, if the knowledge of Avj is readily accessible, one 
may proceed with the estimation of capture probabilities 
in the conventional fashion. If not however, it is sensible to 
assume that this difference of angles is uniformly distributed 
between 0 and 27r, and correspondingly average over Atu, 
leading to the following dehnition of the apsidally-averaged 
capture probability: 

{V) = ^ Vdl^w. (45) 

With this dehnition at hand, the number of input param¬ 
eters is reduced to four (specihcahy, they are the total 
secondary mass (mi + m 2 )/M*, the secondary mass-ratio 
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Figure 8. A series of maps of the apsidally averaged capture probability corresponding to the 2:1 mean motion resonance, for a variety 
of planet-planet and planet-star mass ratios. The color scale marks the averaged capture probability, with orange and blue corresponding 
to high and low values respectively. Additionally, on each plot a {V) = 100%, curve below which capture is certain is shown with a thick 
red line, while a series of contours corresponding to {V) = 75%, 50%, 35%, 25%, 15%, 10% are shown as gold curves. Note that the range 
of the plots increases with total secondary to primary mass. Specifically, emax = 0.08,0.2,0.35 for (mi + m 2 )/M* = 10“^, 10“"^, 10“^ 
respectively. 


and the two orbital eccentricities ei and 62 ) and can 
thus be represented on a series of 2D level contour plots. 

We have calculated the averaged probability across a 
range of parameters for the 2:1, 3:2 and the 4:3 mean motion 
resonances. In particular, we have mapped out {V) on a 
chain of ei vs 62 diagrams, corresponding to total secondary 
masses of (mi +m 2 )/M* = 10 “^, 10 “^, 10 “^ and secondary 
mass-ratios of mil m 2 — 1/9,1,9. The results are shown in 
Figures ( 8 ) - (10). 

On each diagram, the orange-blue color-scale denotes 


the averaged capture probability, which is also marked by 
distinct yellow contour lines. A thickened red line that rep¬ 
resents the conditions for guaranteed capture is also shown 
in each panel. Importantly, within the section of parameter 
space restricted by the red line, any combination of orbital 
eccentricities and differences in apsidal lines will necessar¬ 
ily lead to resonant locking. Note also that the range over 
which ei and 62 are plotted diminishes with decreasing total 
secondary mass. 

Upon examination of Figures ( 8 ) - (10), a few features 


© 0000 RAS, MNRAS 000 , 000-000 


















12 K. Batygin 


3:2 MMR 


mi 

m2 


1/9 


1 


9 



Figure 9. Same as Figure (8), but for the 3:2 mean motion resonance. 


immediately stand out. First, the overall capture probability 
decreases with /c, meaning that capture into the 2:1 mean 
motion resonance is the most generous with respect to the 
degree of pre-encounter orbital excitation. A second read¬ 
ily apparent quality is the weak dependence of {V) on the 
distribution of masses among the two bodies. That is, the 
difference between the values of eccentricities that charac¬ 
terize guaranteed capture is smaller than a factor of two, 
depending on whether the majority of the mass resides in 
the inner or the outer planet. Moreover, the dependence on 
mil m 2 further subsides with increasing k. Third, it is clear 
that the probability of capture most strongly depends on 
the total planetary mass of the system. 

The construction of capture probability maps in prin¬ 
ciple removes the necessity for computationally expensive 


numerical experiments, provided that the assumptions in¬ 
herent to the theory are satisfied for the problem at hand. 
The results also highlight an important attribute of reso¬ 
nant dynamics: while substantial orbital excitation (e ~ 0.1) 
does not inhibit capture for giant planets, rather minute 
(e ~ 0.02) eccentricities may be sufficient to obstruct the 
onset of resonant evolution of Earth-like terrestrial planets. 
We shah return to this point again below. 


3.1.2 Critical Rates of Convergent Evolution 

An implicit assumption inherent to the calculation of cap¬ 
ture probability maps depicted in Figures (8) - (10) is that 
orbital convergence is much slower than angular momentum 
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Figure 10. Same as Figure (8), but for the 4:3 mean motion resonance. 


and energy exchange via resonant interaction. As shown in 
the previous section, it is possible to write down an analytic 
criterion for adiabatic evolution. Accordingly, we shall now 
apply this criterion to the most common mechanisms re¬ 
sponsible for orbital migration (namely interactions of plan¬ 
ets with gaseous disks, scattering of planetesimals, and tidal 
evolution), and derive a series of constraints on physical pa¬ 
rameters inherent to each physical setting. 

• Disk-Driven Migration 

Gravitational interactions between a planet and a 
gaseous disk are well-known to yield a time-irreversible ex¬ 
change of angular momentum (Goldreich & Tremaine 1980). 
At the detailed level, the rate and direction of migration 
are determined by the physical properties of the disk (e.g. 


density and entropy profiles) as well as the planetary mass 
(Paardekooper et al. 2010; Paardekooper, Baruteau, & Kley 
2011; Bitsch & Kley 2011). Depending on the latter, disk- 
driven orbital migration typically falls into one of two^^ dis¬ 
tinct regimes (Ward 1997; Papaloizou & Larwood 2000). 
For planets that are not sufficiently massive to clear out 
their own orbital neighborhoods and thereby carve out gaps 
in their gaseous disks, migration proceeds in the so-called 
“type-I,” or linear regime. We consider this mode of orbital 
transport first. 

A rudimentary way to define a planetary mass below 


For the purposes of this work, we shall neglect the so-called 
“type-III” mode of orbital transport. 
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Figure 11. The adiabatic threshold for two objects undergoing 
type-I migration in an isothermal disk with a Mestel-like surface 
density profile. While this particular figure shows the adiabatic 
limit for a secondary mass-ratio of mi/m 2 = 2/3 = 3/5), 

it should be kept in mind that unlike the expression for rub, 
equation (48) does depend strongly on this quantity. 


which type-I migration applies is to consider the viscous 
gap-opening criterion: 





(46) 


where h is the scale-height of the disk, and a is the Shakura- 
Sunayev viscosity parameter (Armitage 2010). The quoted 
expression is derived by equating resonant and viscous 
torques in the planetary neighborhood. A more sophisti¬ 
cated treatment of the gap-opening process can be found 
in (Crida, Morbidelli, & Masset 2006). 

Under the simplifying assumptions of a power-law sur¬ 
face density profile and an isothermal equation of state, the 
direction of type-I migration is strictly inwards, and occurs 
on a characteristic timescale (Tanaka, Takeuchi, & Ward 
2002 ): 


1 rhV I 

s mTi \a J n’ 


(47) 


where s is a dimensionless constant that depends on the 
structure of the disk. 

Motivated by observations of Andrews et al. (2013), 
here we shall assume that the surface density profile takes 
the form of a classical Mestel disk i.e. S = Soao/a, where 
So is the surface density at some reference semi-major axis 
ao (Mestel 1963). With this choice for S, the have s = 3.8 
(Tanaka, Takeuchi, & Ward 2002). Additionally, for sim¬ 
plicity we shall assume that the disk aspect ratio, /i/a, is 
constant throughout the region of interest. 

From the form of equation (47), it is immediately clear 
that if two planets are undergoing type-I orbital decay, a nec¬ 
essary requirement for convergent migration is m 2 > C^i- 
Assuming that this criterion is satisfied, we evaluate the de¬ 
cay timescales (47) at nominal resonance, such that ai = 
{{k — l)/kY^^ — <^ 2 , and employ equation (44) to obtain 

an expression for the adiabatic threshold: 


3 ^/^ [jrll{mi + 7722 )^ 

(^{{k — ly ky^^ mi — {k — l)?Ti2^ 


(48) 


where Mdisk = Stt Sq uq <22 is the disk mass contained interior 
to the outer planetary orbit^^. This relationship is depicted 
for a particular mass ratio and disk aspect ratio in Figure 

( 11 ). 

It is worth noting that while the resonant dynamics 
themselves tend to only depend on the cumulative plane¬ 
tary mass, the above expression is indeed strongly dependent 
on the mass ratio. The additional dependence on the disk 
surface density profile and separation from the central star 
yields a diverse range of physical parameters which may fall 
either into the strongly adiabatic or strongly non-adiabatic 
regimes (see e.g. Papaloizou & Szuszkiewicz 2005; Cresswell 
& Nelson 2008; Ketchum, Adams, & Bloch 2011). 

Having considered convergent migration of sub-Jovian 
planets, let us now examine the case where the planetary 
mass is sufficient to gravitationally clear out a substantial 
gap around its co-orbital neighborhood. In this situation, 
the planet is shepherded towards the central part of the gap 
where all torques instantaneously cancel. Maintaining this 
configuration, the planet drifts inwards along with the ac¬ 
cretionary flow of the gas (Lin & Papaloizou 1986; Crida & 
Morbidelli 2007; Morbidelli & Crida 2007). Therefore, the 
orbital evolution rate in this so-called ”type-II” regime is 
largely controlled by the global angular momentum trans¬ 
port within the disk^^, and correspondingly operates on the 
viscous timescale: 


2^ _ _2 1 l'h\ ^ 1 
3 ly 3 a \a J n 


(49) 


Here, ly is the kinematic viscosity of the gas (Armitage 2010). 

Equation (49) implies that within the context of the en¬ 
visioned (simplified) picture, two planets that are simulta¬ 
neously undergoing type-H migration will never experience 
convergent resonant encounters. Instead, convergent evolu¬ 
tion is only possible if a less massive object (that migrates 
in the type-I regime) on an external orbit catches up to a 
more massive, slowly migrating planet. For this scenario to 
hold true, the outer planet (m 2 ) must fail the gap-opening 
criterion (46). 

For resonant capture to be successful within the frame¬ 
work of this setup, the system must satisfy two additional 
criteria. Specifically, the outer body must not be so small 
that its type-I orbital decay proceeds at a slower rate than 
type-H migration of the inner body. Simultaneously, the 
outer body must not be massive enough for its migration 
rate to overwhelm the adiabatic limit (43). These two re¬ 
quirements bracket the product of m 2 and Mdisk from below 


In this definition, it is implicitly assumed that a 2 greatly ex¬ 
ceeds the truncation radius of the disk. 

Strictly speaking, this is only true if the planet mass is not 
large enough to disrupt the process of viscous accretion. In other 
words, the dominant portion of the angular momentum budget 
within the planet’s neighborhood must reside in the disk mate¬ 
rial. Quantitatively, this criterion is satisfied when m < 47rEa^, 
where E is the gas surface density immediately outside the gap 
(Baruteau et al. 2014). 
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Figure 12. The adiabatic threshold for two objects, where the in¬ 
ner planet has the planet-star mass ratio identical to Jupiter and 
migrates in the type-II regime, while the outer planet of mass 
m 2 undergoes type-I migration in an isothermal Mestel-like disk. 
The red, blue, and green curves labeled by resonant period ratios 
correspond to the adiabatic criteria and are given by the RHS of 
equation (50). The intersection between the cyan line (given by 
the LHS of equation 50) and the gray curve of slope 2 denotes 
the minimum mass, below which type-II migration proceeds at 
a faster rate than type-I migration, leading to divergent evolu¬ 
tion. The assumed disk parameters are quoted in the figure and 
roughly mimic those of the minimum mass solar nebula. Note that 
for the utilized fiducial parameter choices, Saturn fails the adia¬ 
batic criterion for the 2:1 resonance with Jupiter, but satisfies the 
threshold that corresponds to the 3:2 resonance. However, cap¬ 
ture into the 2:1 resonance is possible, given a diminished surface 
density profile. 

and above: 

Stt k ^ ( m2\( Mdisk \ 

s fc -1 (aj M* j 

A graphic representation of this criterion for a Jupiter-mass 
inner planet is shown in Figure (12). 

In obtaining equation (50), we have made a series of 
approximations. In particular, we used a simplified form of 
the adiabatic threshold as before. Subsequently, we assumed 
that the majority of the planetary mass is contained within 
the inner body. Finally, for the sake of a less cumbersome 
expression, we assumed that convergent evolution is domi¬ 
nated by type-I migration of the outer body in the second 
inequality. These simplifications generally hold as long as 
the bracketed quantity is not too close to either limit. 

For typical disk parameters (e.g. h/a ~ 0.05; a ~ 
0.001), the extrema in equation (50) differ by a little more 
than an order of magnitude. However, it should be noted 
that without invoking unreasonable quantities, it is possible 
to make the range of planetary masses for which adiabatic 
capture may take place, exceptionally small. 

• Planetesimal-Driven Migration 

In the post-nebular stage of planetary system evolu¬ 
tion, torques derived from interactions with the gas are non¬ 
existent. However, large-scale migration can still occur as a 
consequence of asymmetric scattering of planetesimals by a 
more massive object (Fernandez & Ip 1984; Malhotra 1995). 



/i 


Figure 13. The adiabatic threshold for planetesimal driven mi¬ 
gration. As in previous figures, the adiabatic criteria are shown 
with red, blue, and green lines. However, these curves lie well be¬ 
low the mass-ratio criterion (shown with a cyan line) which must 
be exceeded for migration to be self-sustaining. This means that 
planetary embryos that migrate through long-lived planetesimal 
disks generally do not suffer adiabatic resonant encounters. 

Within the context of the early dynamical evolution of the 
outer Solar System, this process may have played a crucial 
role, as scattering-facilitated divergent evolution of Jupiter 
and Saturn has been invoked as a means of igniting a tran¬ 
sient dynamical instability^^ (Tsiganis et al. 2005; Batygin 
& Brown 2010). 

A prominent setting where planetesimal-driven migra¬ 
tion can lead to substantial orbital changes is a gas-depleted 
planetesimal disk with embedded planetary embryos. In 
light of this, it is worthwhile to examine if resonant lock¬ 
ing among such embryos can occur. A fiducial timescale on 
which planetesimal-driven migration operates is (Ida et al. 
2000; Kirsh et al. 2009) 


where E now refers to planetesimal surface density, for which 
we have once again assumed a Mestel-like radial profile. A 
detailed study of planetesimal-driven migration was recently 
carried out by Minton & Levison (2014), who identified five 
criteria that must be satisfied for this mode of orbital trans¬ 
port to be self-sustained. Here, we shall not dwell extensively 
on these requirements, but simply assume that they are sat¬ 
isfied for at least one of the secondary objects of interest. 

The direction of planetesimal-driven migration can be 
inward or outward. To this end, Kirsh et al. (2009) find a 
preference for orbital decay over growth, however Minton & 
Levison (2014) find no statistically significant predisposition 
for either direction. From equation (51), it is clear that the 
inner orbit changes on a shorter timescale. Therefore, if both 
orbits decay at rates given by equation (51), the evolution 
is bound to be divergent. Reversing this argument, one may 

Among the many features whose origins can be attributed to 
this instability are the so-called period of late heavy bombard¬ 
ment (Gomes et al. 2005; Levison et al. 2011), Jupiter’s Trojan 
asteroids (Morbidelli et al. 2005), irregular satellite populations of 
the giant planets (Nesvorny, Vokrouhlicky, & Morbidelli 2007), 
and the dynamical structure of the Kuiper belt (Levison et al. 
2008; Batygin, Brown, & Fraser 2011). 
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envisage that convergent evolution is possible if both orbits 
grow. However, in this case the inner planet would be mi¬ 
grating in the wake of the outer planet, through a pre-excited 
disk. Without invoking very rapid dissipation that would dy¬ 
namically cool the disk, such a scenario would likely violate 
the disk eccentricity criterion (see Minton & Levison 2014), 
halting the inner orbit’s migration. Consequently, it would 
appear that given nominal parameters, convergent evolution 
is possible either if the inner orbit evolves outwards while 
the outer evolves inwards, or if one of the orbits remains 
stationary while being approached by the other. 

In principle, each of these situations will be character¬ 
ized by its own, subtly different adiabatic threshold. How¬ 
ever, in interest of succinctness we shall only consider the 
case where both orbits actively approach one-another, not¬ 
ing that the corresponding expressions for the case where 
one of the planets remains stationary are quantitatively sim¬ 
ilar, to within a factor of ~ 2. Combining expressions (51) 
and (43), we have: 


MdiA < 3"/=^ 

M* ~ TT fc-l + (fc(fc-l)2)l/3 
( /reV (mi + m2) ^ 

"" V ) ■ 


(52) 


It is interesting to note that while expression (52) pro¬ 
vides an upper bound on the disk mass, the mass ratio 
criterion for operation of planetesimal-driven migration re¬ 
quires the mass of the migrating body to not exceed its 
encounter mass (equivalently, isolation mass; Lissauer 1987) 
by more than a factor of ~ 3 (Kirsh et al. 2009), thereby 
providing a lower bound. Given that the multiplicative fac¬ 
tor of (mi + m 2 )/Mi, is of order unity for all resonances 
of interest and that the stopping mass is approximately 
(m/M*)stop — 12 (Mdisk/M*)^/^, we can formulate the fol¬ 
lowing rough criterion for scattering-facilitated adiabatic 
resonant encounters: 


Mdisk ^ / mi + m 2 ^ ^ / Mdisk \ ^ 

M* ~ V M* J ~ V / 


(53) 


Despite its crudeness, the above expression is informa¬ 
tive, and is shown in Figure (13). In particular, it dictates 
that in order to simultaneously satisfy the adiabatic crite¬ 
rion and the mass ratio criterion, the planetesimal disk mass 
must exceed Mdisk/M* > 0.04. Such a mass ratio is charac¬ 
teristic of gaseous protoplanet ary disks, and thus exceeds the 
mass ratio inherent to the solid component of the disk by 
approximately two orders of magnitude. Therefore, equation 
(53) implies that without appealing to special configurations 
(such as the primordial multi-resonant state of the outer So¬ 
lar System), planetary migration facilitated by scattering of 
planetesimal cannot lead to convergent adiabatic resonant 
encounters. 

• Tidal Evolution 

A final mode of convergent migration that we shall con¬ 
sider here is that facilitated by tidal dissipation. This mecha¬ 
nism is particularly relevant to planetary satellites, as orbital 
changes induced by tides raised on the host planet are con¬ 
sidered to be the dominant driving mechanism^® responsible 


See however, Peale & Lee (2002) and Canup & Ward (2002) for 





Figure 14. The adiabatic threshold corresponding to tidal evolu¬ 
tion of an equal-mass pair of satellites. For illustrative purposes, 
in this figure we considered the inner object’s orbit to approxi¬ 
mately coincide with the Roche radius, yielding a maximal mi¬ 
gration rate. Even in this extreme case, an reasonable dissipation 
efficiency will yield adiabatic evolution. 


for resonant pairs of satellites in the Solar System (Goldre- 
ich 1965; Peale 1976, 1986; Henrard 1983). Correspondingly, 
we shall frame the following discussion in the satellite-planet 
context, keeping in mind that similar arguments can apply 
to close-in planets and their host stars (see e.g. Adams & 
Bloch 2015). 

Following the works of Allan (1969) and Dermott, Mal- 
hotra, & Murray (1988), we shall assume that changes in the 
semi-major axis induced by satellite tides can be neglected 
in favor of their planetary counterparts. In this case, the 
characteristic migration timescale for an orbit residing be¬ 
yond the synchronous radius is given by (Murray & Dermott 
1999): 


1 Q ^ mV - 

3 /c2 m \Rp J n ’ 


(54) 


where Q is the specific dissipation function, k 2 is the Love 
number, and Rp is the physical radius of the planet. 

Contrary to the case of type-I migration discussed 
above, the expression (54) implies that the outer satellite 
mass cannot exceed the inner satellite mass by more than 
m 2 < mi(k/k — 1)^^/^ for evolution to remain convergent. 
Substitution of the tidal timescale into equation (44) yields 
the following expression for the adiabatic threshold: 


Q ( (^2 
k2 \Rp 




M, 


2 \fielimi+m2) 
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f - {k - 

V (fc - l)4i/9Mp 


(55) 


To evaluate the physical meaning of equation (55), it 
is instructive to consider a limiting parameter regime that 
maximizes the migration rate. Accordingly, let us envision 
that the total secondary mass is concentrated entirely in the 
inner satellite which orbits immediately outside the Roche 
limit (i.e. 02 — (5/2) Rp), and the planet is a homogenous 
fluid body with a Love number of M = 3/2. Setting the 


an alternative, circumplanetary disk driven view on the assembly 
of the Laplace resonance among the Galilean satellites. 

© 0000 RAS, MNRAS 000 , 000-000 












MMR Capture Theory 17 


speciflc dissipation function to a minimal^^ value of Q ~ 10, 
we can obtain a critical mass ratio (mi/Mp)crit below which 
the adiabatic limit is broken. Upon evaluation, we obtain 
(mi/Mp)crit ~ and 10“^ for the 2:1, 3:2, and 

4:3 resonances respectively. 

While these values already lie well below the typi¬ 
cal range of natural^® satellite - planet mass ratios, it is 
important to keep in mind that we have invoked unrea¬ 
sonably rapid tidal evolution (as an example, for Jupiter, 
k 2 /Q ~ 10“^ whereas we have used /C 2 /Q ~ 3/20; Gol- 
dreich & Soter 1966; Lainey et al. 2009) in order to de¬ 
rive them. In other words, the obtained quantities are gross 
overestimates of the practically relevant critical mass ratios. 
Recall also that the rate of tidal evolution diminishes with 
(m/Mp) meaning that the critical mass ratio may also lie in 
a regime where tt greatly exceeds the age of the system, ren¬ 
dering this quantity meaningless. Therefore, for all practical 
purposes it is likely safe to assume that tidally-facilitated 
resonant encounters will almost always lie in the adiabatic 
regime. For completeness. Figure (14) shows the adiabatic 
threshold for an equal-mass pair of satellites. 

The calculations presented above yield a series of prac¬ 
tically useful criteria that inform whether or not convergent 
migration within a given physical setting adheres to the adi¬ 
abatic limit. In interpreting this discussion however, it is 
important to keep in mind that for the sake of definitive¬ 
ness, we have limited our calculations to simplified formulae 
that describe nominal dissipative evolution rates of planets. 
In reality, the detailed physics of disk-planet interactions, 
planetesimal scattering and tidal evolution can be rather 
complex (see Kley & Nelson 2012; Baruteau et al. 2014; 
Minton & Levison 2014; Efroimsky & Williams 2009 for in- 
depth discussions). Indeed, substantial deviations from the 
quoted prescriptions can occur within a more comprehensive 
treatment. As a result, the derived criteria should be viewed 
as approximate, rather than conclusive. 

3.2 Specific Applications 

3.2.1 Resonant Capture of Jupiter and Saturn in the 
Protosolar Nebula 

As already mentioned above, type-II migration generally 
causes planets to spiral in towards their host stars. In light of 
this, the following question naturally arises: why are Jupiter 
and Saturn where they are today? As an answer to this ques¬ 
tion, it was shown by Masset & Snellgrove (2001) that the 
simultaneous dynamics of Jupiter and Saturn, when sub¬ 
merged in the protosolar nebula, alter the isolated migra¬ 
tion picture qualitatively. In particular, simulations show 
that Jupiter and Saturn exhibit convergent orbital evolu¬ 
tion, and eventually get captured in resonance. With a res- 

The smallest possible value that Q can take on is unity. Phys¬ 
ically, this would correspond to complete dissipation of all energy 
stored within a single tidal cycle. However, for the weak-friction 
theory (within the context of which equation (54) is obtained) to 
supply, Q must greatly exceed unity (Hut 1981). 

Note that the assumptions of low eccentricity and inclination 
inherent to the formulation of the Hamiltonian (1) render our 
treatment largely inapplicable to irregular satellites of the Solar 
System. 


onant lock established, the planets carve out a mutual gap, 
which (owing to a particular planet-planet mass ratio) al¬ 
ters the torque balance of the system and facilitates joint 
outward migration of the resonant pair (Morbidelli & Crida 
2007). 

Although resonant reversal of orbital decay was initially 
proposed to account for the retention of Jupiter and Saturn 
at large orbital radii, such a sequence of events smoothly 
connects the nebular stage of Solar System evolution to the 
compact initial conditions of the Nice model (Morbidelli et 
al. 2007; Batygin & Brown 2010; Nesvorny & Morbidelli 
2012 ), and provides natural explanations for the compara¬ 
tively small mass of Mars (Walsh et al. 2011), the composi¬ 
tion of the Asteroid belt (O’Brien et al. 2014), and the non¬ 
existence of otherwise common close-in Super-Earths inside 
of Mercury’s orbit (?). Moreover, it has been proposed that 
a similar mechanism is ubiquitously responsible for halting 
large-scale migration of giant exoplanets that reside at wide 
orbital separations (Morbidelli 2013). 

Eollowing the pioneering study of Masset & Snellgrove 
(2001), early dynamical evolution of Jupiter and Saturn has 
been studied by Morbidelli & Crida (2007); Pierens & Nelson 
(2008); D’Angelo & Marzari (2012); Pierens et al. (2014) us¬ 
ing hydrodynamical simulations. In agreement with the ini¬ 
tial experiments of Masset & Snellgrove (2001), most stud¬ 
ies find that independent of initial conditions, the common 
outcome is passage through the 2:1 resonance and capture 
into the 3:2 resonance. The prevalent speculation (and in¬ 
deed this turns out to be correct) within the aforementioned 
works is that failure of the 2:1 resonance to capture Jupiter 
and Saturn arises from non-adiabatic evolution. Given the 
well-defined nature of the problem, it is instructive to ex¬ 
amine how the results of numerical simulations ht into the 
framework of the developed analytic theory. 

As a first step, let us imagine that the adiabatic crite¬ 
rion is satisfied. In this regime, we can use equation (22) to 
inform the conditions under which resonant capture is not 
guaranteed. Specific curves corresponding to the 2:1 Jupiter- 
Saturn resonance are shown in Eigure (7). Evidently even in 
the conservative case, pre-encounter eccentricities in excess 
of e > 0.1 are required to reduce the adiabatic capture prob¬ 
ability below ~ 75%. Gonventional theory (e.g. Goldreich & 
Tremaine 1980; Ward 1997) suggests that dissipative evo¬ 
lution of planets within axisymmetric disks tends to damp 
orbital eccentricities and inclinations^^, meaning that cap¬ 
ture into the 2:1 resonance is highly likely if the adiabatic 
limit is satisfied. Recall also that the 2:1 resonance allows 
for guaranteed capture at higher eccentricity than the 3:2 
resonance, meaning that if the system fails to capture into 
the 2:1 resonance due to being too dynamically pre-excited, 
it will also likely fail to capture into the 3:2 resonance. 

Let us now evaluate the critical rates of convergent evo¬ 
lution for Jupiter and Saturn in the protosolar nebula, as 
dictated by equation (50). Adopting the typically quoted as¬ 
pect ratio of h/a = 0.05, a viscosity parameter of o: = 10~^, 
a surface density of So = 2000 g/cm^ at ao = 1 AU, and 
an encounter semi-major axis of a 2 — 2.5 AU (the envi¬ 
sioned “tacking” radius), the minimum mass of the outer 


See Goldreich & Sari (2003); Tsang (2014) for an alternative 
view. 
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planet above which the orbits will approach each-other is 
m 2 > 3 M 0 . This is in excellent agreement with the simula¬ 
tion results of Pierens & Nelson (2008) who find the mini¬ 
mum mass to be m 2 ^ 3 . 5 M 0 . 

Using the same expression we find the threshold mass, 
above which the adiabatic criterion is violated for the 2:1 
resonance, is m 2 ^ 6 OM 0 , or about 2/3 of Saturn’s mass. 
For the 3:2 resonance, the critical value increases to m 2 > 
I 4 OM 0 , which comfortably exceeds Saturn’s mass. Again, 
these estimates are in perfect alignment with the numerical 
experiments of Pierens & Nelson (2008) who observe capture 
into the 2:1 resonance for m 2 = 30 — 4 OM 0 and capture into 
the 3:2 resonance for m 2 = 80 — IOOM 0 . Thus, our analyti¬ 
cal theory robustly conforms to the results of numerical ex¬ 
periments and analytically demonstrates that the recurrent 
capture of Saturn into a 3:2 rather than a 2:1 mean mo¬ 
tion resonance is indeed an outcome of the violation of the 
adiabatic threshold for the 2:1, but not the 3:2 resonance. 

Equation (50) further implies that capture into the 2:1 
resonance is indeed possible, provided a reduced disk surface 
density. As an example, in addition to results correspond¬ 
ing to fiducial parameters quoted above. Figure (12) also 
shows a supplementary (dashed) curve that represents a disk 
with Eo = 900 g/cm^. The fact that capture is permitted in 
this regime is fully consistent with the simulation results of 
Pierens et al. (2014), who numerically obtain capture into 
the 2:1 resonance for sub-nominal surface densities. 

Naturally, the analytic model developed herein does not 
take into account the full richness of the possibilities that can 
occur in real protoplanetary disks. For example, simulations 
show that after Saturn passes through the 2:1 resonance, 
it directly invades the disk gap opened by Jupiter. As a 
result, the diminished local surface density of the gas further 
reduces its migration rate, allowing for an even more robust 
satisfaction of the adiabatic criterion (Morbidelli & Crida 
2007). A similar reduction in Saturn’s migration rate can in 
principle occur if Saturn successfully opens its own gap, such 
that its migration rate is still faster than the strict type-II 
rate given by equation (49), but only slightly so. As is made 
clear by equation (46), this can be achieved by reducing the 
aspect ratio or the viscosity of the disk. 

3.2.2 Paucity of Resonances Among Close-In Sub-Jovian 
Planets 

Having examined the well-studied problem of nebular evo¬ 
lution of Jupiter and Saturn, let us now consider a some¬ 
what less transparent phenomenon. A pivotal realization 
that has emerged as a consequence of the recent identifica¬ 
tion of thousands of extrasolar planets by the Kepler tran¬ 
sit survey, is that the dominant mode of planet formation 
within the Galaxy is one that produces planets substantially 
smaller than Jupiter and Saturn with orbital periods below 
~ 100 days (Youdin 2011; Howard et al. 2012). These close- 
in sub-Jovian planets are estimated to orbit approximately 
half of all Sun-like stars, and often comprise tightly packed 
multi-planet systems (Mayor et al. 2011; Batalha et al. 2013; 
Lissauer et al. 2014). Accordingly, understanding the archi¬ 
tectural origins of such systems holds the key to identifying 
the dominant processes at play during the epoch of planet 
formation. 

The orbital distribution of close-in sub-Jovian planets is 


for the most part, devoid of orbital resonances. While there 
exist substantial enhancements in the periods ratios imme¬ 
diately outside of the 2:1 and 3:2 commensurabilities (and 
such systems can indeed be attributed to resonant evolution 
in presence of dissipative forces; Lithwick & Wu 2012; Delisle 
et al. 2012; Batygin & Morbidelli 2013a), these planet-pairs 
comprise a minority within the observational sample (Fab- 
rycky et al. 2014). This fact is puzzling when viewed within 
the context of the theoretical expectation that type-I mi¬ 
gration should typically generate crowded resonant chains 
in protoplanetary disks (Cresswell & Nelson 2008). 

One way to prevent resonant capture is to invoke a suf¬ 
ficiently massive disk such that the adiabatic criterion inher¬ 
ent to type-I migration (equation 48) is not satisfied (Quillen 
2006; Mustill & Wyatt 2011). Although this will indeed pro¬ 
vide an effective avenue towards halting the production of 
resonant chains of planets, this process alone is unlikely to 
explain the observed exoplanet architectures in a satisfac¬ 
tory manner for two reasons. First, the disk mass itself de¬ 
creases in time, generally on a timescale that is much longer 
than typical type-I migration timescales (Haisch, Lada, & 
Lada 2001; Herczeg & Hillenbrand 2008). This means that 
even if a given disk starts out in a state where resonant 
capture is hampered by overly rapid migration, such a disk 
will eventually evolve to a state where the resonant condi¬ 
tion (48) is satished. Second, it has been shown that even in 
massive disks, type-I migration is not ubiquitously fast, as 
suggested by conventional linear calculations (Paardekooper 
et al. 2010; Paardekooper, Baruteau, & Kley 2011; Bitsch & 
Kley 2011; Bitsch, Boley, & Kley 2013). The main deviation 
from linear theory arises from the fact that entropy gradients 
in radiative disks alter the disk-planet interaction in such a 
way as to create zones where type-I migration greatly slows 
down, reverses, or ceases all together, allowing for conver¬ 
gent congregation of low-mass planets (McNeil, Duncan, & 
Levison 2005; Papaloizou & Szuszkiewicz 2005; Matsumoto, 
Nagasawa, & Ida 2012). 

An alternative means of generating a largely non¬ 
resonant orbital distribution is to invoke processes that tend 
to destabilize mean-motion commensurabilities. To this end, 
Adams, Laughlin, & Bloch (2008), followed by Ketchum, 
Adams, & Bloch (2011) and Rein & Papaloizou (2009) pro¬ 
posed that stochastic perturbations due to turbulent forc¬ 
ing within protoplanetary disks may inhibit long-lived res¬ 
onances (see also Laughlin, Steinacker, & Adams 2004 for 
a related discussion). Another possibility, recently put for¬ 
ward by Goldreich & Schlichting (2014) is that the com¬ 
bined effects of semi-major axis and eccentricity damp¬ 
ing may render resonant configurations unstable on long 
timescales. While both of these proposals are in principle 
reasonable, they ultimately rely on processes that are poorly 
constrained. Specifically, these mechanisms inherently de¬ 
pend on quantities such as the duty-cycle of turbulent eddies 
and the relative damping rates of orbital eccentricities and 
semi-major axes, both of which remain scantily understood 
and constitute active fields of research (Bai & Stone 2013; 
Bitsch et al. 2015). As a result, in the spirit of Occam’s ra¬ 
zor, it is tempting to inquire if a simpler mechanism that 
can preclude resonant capture exists. 

Keeping in mind the arguments presented above, let us 
only consider the adiabatic regime. As we already saw in the 
beginning of this section, capture is guaranteed only below 
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Figure 15. Distributions of critical orbital eccentricities, above which resonant capture of confirmed planets pairs of the Kepler survey is 
not guaranteed. Because capture probability decreases rapidly with pre-encounter orbital excitation (see Figure 3B), the quoted average 
eccentricities are approximately indicative of the true eccentricities at which the outcomes of adiabatic resonant encounters switch from 
preferential capture to preferential transition across the resonance. Results presented in panels A, B, C correspond to 2:1, 3:2, 4:3 mean 
motion resonances respectively. 


some critical extent of pre-encounter orbital excitation, and 
capture probability decreases rapidly above this value (Fig¬ 
ure 3B). How eccentric must have close-in planet pair been 
to preferentially skip over resonances without locking? To 
obtain a rough answer to this question, consider the follow¬ 
ing calculation. 

Suppose that by default, planets in protoplanetary disks 
migrate not on exactly circular orbits but with some small 
eccentricities, e. Further, let us further speculate that e will 
be similar for two planets that are about to encounter a 
mean motion resonance, as it is set by external factors rather 
than direct planet-planet interactions. Provided that plane¬ 
tary and stellar masses are known, the critical value of pre¬ 
encounter eccentricity, dcrit, can be computed by means of 
equation (22), setting Atu = tt. We have performed this 
calculation for all confirmed planet-pairs detected by the 
Kepler transit survey. 

Figure (15) shows the distributions of Ccrit correspond¬ 
ing to the 2:1, 3:2 and 4:3 resonances for Kepler plan¬ 
ets. Clearly, the pre-encounter orbital eccentricities required 
to render capture uncertain are ubiquitously low: Ccrit — 
0.01 — 0.03. More importantly however, our simple estimate 
for dcrit is astonishingly close to the observed rms eccentrici¬ 
ties of the Kepler sample. Specifically, from their analysis of 
transit timing variations, Wu & Lithwick (2013) and Hadden 
& Lithwick (2014) find that three quarters of the considered 
sub-sample of planets conforms to an eccentricity distribu¬ 
tion with an rms free eccentricity of e ~ 0.02, while the re¬ 
maining quarter is characterized by substantially higher val¬ 
ues. Although caution must be exercised in interpreting the 
present eccentricities of sub-Jovian planets as having been 
inherited from their natal disks, the conspicuous similarity 
of the theoretically required values to the corresponding ob¬ 
servations is suggestive, and opens a previously unexplored 
avenue for explaining the origins of the orbital distribution. 


4 DISCUSSION 

Both theoretical and observational lines of inquiry suggest 
that mean-motion resonances play a central role in the for¬ 
mation and long-term evolution of planetary systems (see 


e.g. Morbidelli et al. 2007; Rivera et al. 2010; Deck et 
al. 2012; Goldreich & Schlichting 2014 and the reference 
therein). While this fact has been well recognized as a con¬ 
sequence of countless numerical experiments (Quillen 2006; 
Ketchum, Adams, & Bloch 2011; Ogihara & Kobayashi 
2013), until now a comprehensive qualitative understand¬ 
ing of the conditions for capture of planets into resonances 
had remained elusive. In this work, we have taken steps to¬ 
wards formulating a model for capture into resonance in a 
generic, analytical way. Our development provides an under¬ 
lying framework within which the outcomes of numerical ex¬ 
periments can be interpreted. Moreover, we provide a series 
of simple criteria which can be used to inform the outcome 
of dynamical simulations. 

As a practical recipe for theoretical analysis of resonant 
encounters, we propose the following order of calculation: 

• Does the convergent evolution rate of the system satisfy 
the adiabatic condition? For the specific cases of orbital mi¬ 
gration driven by interactions with a gaseous disk, scattering 
of planetesimals, or tidal evolution, expressions (48) or (50), 
(53), and (55) may be used respectively. Alternatively, the 
generalized adiabatic criterion (43), or its simplified form 
(44) can be employed. If the criterion is violated, capture 
will not take place. 

• If the resonant encounter lies in the adiabatic regime, is 
capture certain? This can be assessed using expression (22). 
If knowledge of the differences in the apsidal lines, Atu, is 
available, the guaranteed capture equation can be used di¬ 
rectly. If Avj is not known, it is sensible to assume Attj = tt, 
as this yields the most conservative estimate. 

• If adiabatic capture is not guaranteed, what is the capture 
probability? If Avj is known, the initial condition far away 
from resonance can be related to the location of the unsta¬ 
ble fixed point of the resonance at separatrix-crossing using 
expression (27). The corresponding probability can then be 
calculated using equation (28). If Avj is unspecified, the av¬ 
eraged capture probability given by equation (45) can be 
calculated. 

In this work, we have utilized the developed formalism 
to consider two specific applications. First, we analyzed the 
numerically well modeled evolution of Jupiter and Saturn’s 
orbits in the primordial Solar nebula (Pierens et al. 2014 and 
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the references therein). Within the context of this problem, 
our theoretical arguments confirmed the previously insin¬ 
uated notion (e.g. Morbidelli & Crida 2007) that the ten¬ 
dency of Jupiter and Saturn to lock into the 3:2 rather than 
2:1 resonance is a consequence of Saturn’s rapid migration 
and the associated violation of the adiabatic criterion. We 
subsequently considered the origins of the dominantly non¬ 
resonant orbital distribution of close-in sub-Jovian multi¬ 
planet systems discovered by Kepler (Batalha et al. 2013; 
Fabrycky et al. 2014). Specifically, we showed that if low- 
mass planets generally reside on slightly eccentric orbits (i.e. 
e > 0.02) when submerged in their natal disks, the chances 
of resonant capture are greatly diminished. 

Although the theoretically derived threshold eccentric¬ 
ities (below which resonant capture is hindered) are al¬ 
most identical to the present observationally inferred values 
within the Kepler sample (Wu & Lithwick 2013), we have 
not specified the physical origin of these subtle deviations 
from circular orbits. In principle, there exists a large num¬ 
ber of dynamical mechanisms that may lead to such excita¬ 
tion. Among them are orbital excitations by turbulent forc¬ 
ing (Nelson & Papaloizou 2004; Adams & Bloch 2009) and 
interactions with distant massive planets which may them¬ 
selves experience violent dynamical instabilities within pro¬ 
toplanetary nebulae (Lega, Morbidelli, & Nesvorny 2013). 

An arguably more intriguing idea is that steady-state 
non-zero eccentricities may stem directly from the interac¬ 
tions of planets with their natal disks^°. A body of recent 
literature has shown that contrary to the conventional the¬ 
oretical simplification of perfectly circular protoplanetary 
nebulae, real disks exhibit substantial deviations away from 
axial symmetry and are indeed believed to be globally lop¬ 
sided (Brown et al. 2009; Tang et al. 2012; Casassus et 
al. 2013; Fukagawa et al. 2013; Isella et al. 2013; van der 
Marel et al. 2013; Perez et al. 2014; Bruderer et al. 2014). 
Although observations taken in the dust continuum wave¬ 
length band severely overrepresent the corresponding non- 
axisymmetric over-densities of gas (which holds the vast ma¬ 
jority of the disk mass), the gas disk eccentricities^^ required 
to explain the ALMA observations (Mittal & Chiang 2015) 
may very well be sufficient to also perturb low-mass planets 
onto slightly non-circular orbits. While in this work, we have 
taken initial steps towards exploring this idea, undoubtably 
much additional effort is required to quantitatively assess 
the viability of this hypothesis. 

Our analytical study of resonant capture complements 
a number of recent developments that employ the same inte- 
grable formalism for resonant dynamics. In a closely related 


Note that the timescale of eccentricity modulation due to 
planet-disk interaction is almost certainly much shorter than 
that corresponding to semi-major axis evolution (Goldreich & 
Tremaine 1980; Lee &: Peale 2002). 

Possible mechanisms that may be responsible for maintenance 
of global lopsided modes of disks include disk self-gravity (Dury 
et al. 2008; Mittal & Chiang 2015; see also Touma, Tremaine, 
& Kazandjian 2009 and Batygin 2012 for a related discussion) 
and hydrodynamic forces (Larwood et al. 1996; Xiang-Gruess & 
Papaloizou 2014). Moreover, excitation of disk eccentricities may 
arise from external perturbations by passing or bound stars in star 
formation environments that are well-known to exhibit enhanced 
stellar multiplicity (Duchene & Kraus 2013). 


study (Batygin & Morbidelli 2013b), we showed how reso¬ 
nant dynamics can be represented in an intuitive geometric 
way, and demonstrated that divergent resonant encounters 
(where capture necessarily fails) leave the system in a per¬ 
sistent apsidally anti-aligned state. Using the same formal¬ 
ism, Deck, Payne, & Holman (2013) considered the onset of 
chaotic motion in the unrestricted elliptic three-body prob¬ 
lem, and showed that the first order resonance overlap crite¬ 
rion is roughly independent of the planet-planet mass ratio. 
This is perfectly congruent with our finding that the process 
of resonant capture and the associated adiabatic threshold 
only exhibit strong dependence on the ratio of the cumula¬ 
tive secondary mass to the primary mass. 

There is a number of ways in which our theory can be 
expanded upon. Accordingly, we wish to conclude the paper 
with a list of possible directions for future development. 

• To obtain a better relationship between the specified or¬ 
bital conditions far away from resonance and those at the 
resonant encounter, it may be fruitful to characterize the 
pre-encounter secular evolution. This may better inform the 
true value of the action at the time of the appearance of 
the separatrix. 

• A more complete version of the model may incorporate the 
effects of external eccentricity damping or excitation (e.g. 
Lee & Peale 2002; Goldreich & Schlichting 2014). Direct 
modulation of the eccentricities will lead to corresponding 
changes in the value of the action and may cause J to shift 
from a value smaller than J7c to a value that exceeds it, 
or vice-verse. It is possible that the inclusion of this effect 
will not affect the formulation of capture probabilities be¬ 
cause one can envision defining time-dependent canonical 
coordinates, in which is conserved by construction and 
the evolution rate of S is correspondingly accelerated or di¬ 
minished (Henrard 1993). However, this assertion deserves 
to be tested explicitly. 

• In our analysis of probabilistic capture, we have adopted 
the conventional approach of considering the crossing of a 
well-defined separatrix, and thus neglected the effects of 
chaos. In fact, retention of higher order terms in eccentricity 
and inclination will increase the number of degrees of free¬ 
dom, rendering the Hamiltonian non-integrable. In a system 
with multiple degrees of freedom, the vicinity of the sepa¬ 
ratrix may be engulfed in a chaotic layer, which may alter 
the outcome of resonant encounters (attempts at character¬ 
ization of crossing of a stochastic layer have been previously 
made by Henrard & Morbidelli 1993). 

• At present, resonant capture theory does not account for 
external stochastic forces. Such effects have however been 
shown to compromise the locking and longevity of reso¬ 
nances (Rein & Papaloizou 2009; Paardekooper, Rein, & 
Kley 2013). A more complete theoretical framework for this 
reduction in capture probability may perhaps be constructed 
with the aid of stochastic calculus (see e.g. Adams, Laugh- 
lin, & Bloch 2008; Batygin, Morbidelli, & Holman 2015 for 
a related discussion). 

• A somewhat more approximate (but similar in spirit to 
what has been done here) integrable model for second-order 
resonances has recently been presented by Delisle, Laskar, & 
Correia (2014). Accordingly, the existing capture theory for 
second order resonances in the restricted problem (Henrard 
1982; Borderies & Goldreich 1984) can be adopted for the 
elliptic problem as well. 
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Cumulatively, it is clear that the extent of theoretical 
expansion that can potentially be undertaken is substantial. 
However, developments such as those proposed above will 
surely contribute to the construction of a comprehensive 
model for planetary system formation and evolution, and 
are thus well-motivated. 
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